Topological Insulator and the ^-vacuum in a system without boundaries 
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In this paper we address two questions concerning the eflective action of a topological insulator 
in one and three dimensional space without boundaries, such as a torus. The first is whether a 
uniform 6'-term with = tt is generated for a strong topological insulator. The second is whether 
such a term has observable consequences in the bulk. The answers to both questions are positive, 
but the observability in three dimension vanishes for infinite system size. 



I. INTRODUCTION 

The topological insulators are characterized by insu- 
lating band structures with nontrivial topology, which 
cannot be smoothly deformed back to an atomic insula- 
tor, sometimes under certain discrete symmetryii"— The 
most well-known example is the integer quantum Hall 
effect (IQHE) in two dimensions (2D). Here despite the 
fact that the bulk of the system is gapped, the system 
possesses n gapless chiral edge states where n is the num- 
ber of occupied Landau levels. In addition, if we perturb 
the system using local electric fields, there will be local 
transverse current in the system. This effect is best cap- 
tured by the Chcrn-Simons term J d^xdte^'''''A^9,yAA in 
the effective theory. 

In three dimensions (3D), under time- reversal (TR) 
symmetry there is a Z2 classification for band 
insulators;^ which distinguishes the usual insulator from 
the topological insulator. In contrast to the IQHE states, 
the two classes of insulators are often distinguished only 
on the edge: in the bulk both insulators do not have a 
current response to the applied field. However, there is 
an odd number of spin-filtered gapless edge states, i.e., 
helical Dirac cones, on the surface of the topological insu- 
lator. Another physical effect of the topological insulator 
can be observed when the TR symmetry is locally broken 
on the edge: the Dirac cones would be gapped and result 
in an ^-integer quantum Hall effect on the edgei^ This ef- 
fect also leads to the bulk electromagneto-polarizabilityj^ 

The effect mentioned above can also be understood 
using the effective topological field theoryi^ This theory 
postulates that the band topological insulator can be de- 
scribed by a bulk " 0-term" , 
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with 6 ~ TT. The integrand in eq. ([T]) is a total deriva- 
tive so the equation of motion remains unaltered in the 
bulk. The edge quantum Hall effect can be derived with 
a smooth change of 9 from tt to in space. Finally, the 
gapless spin-filtered edge states, are understood as a limit 
when one preserves the TR symmetry on the boundary. 

There is, however, some confusion as to whether the 
bulk 0-term, given by eq. ^ is indeed required to de- 
scribe the topological insulator. Specifically, the "V6- 



term" 
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which differs from eq. ([T]) by a total derivative, is what 
actually has been derived in the previous theories and can 
produce all the effects mentioned above as well. There 
is an important conceptual distinction between the two 
theories: with the former theory, it looks as if the bulk 
of the topological insulator has some special property, 
which manifests itself either on the edges or when the 
Hamiltonian is varied; with the latter theory, the bulk 
of the topological insulator is not intrinsically different 
from an ordinary insulator, and all the possible physical 
consequences only occur at the edge between it and a 
normal insulator, or with a variation of the Hamiltonian. 

In this paper, first we show that there are physical 
consequences which distinguish between the two theo- 
ries, and then we show that it is the former theory, Cg, 
which describes the topological insulator. The paper is 
organized as follows: in the following section, we study 
the effective theory on a closed manifold and show that 
the two theories given by eq. ([T|) and eq. ([2]) behave dif- 
ferently. In section HI, we start from the fermionic band 
structure and derive the effective theory Ce, without the 
ambiguity of a total derivative. 



II. THE PHYSICAL CONSEQUENCE OF THE 9 
TERM 



Without a boundary, C^/g will have no effect if 9 is 
uniform. Therefore, the distinction between the two afor- 
mentioned theories will be evident if there is any physical 
consequence of a O-tevm with uniform 9. In the following 
we shall discuss the effect of an uniform O-tevm. in one 
and three dimensions, with various topologies. 

9-term in one spatial dimension (ID): 

Here we shall follow the approach of the 0-vacuum, 
where the 9 term gives a prescription to form gauge- 
invariant states. Please see appendix A for a derivation 
directly using the path integral. 

For completeness we first show that the 0- vacuum de- 
scription is equivalent to a path integral with Cg- The 
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Let us first take the = gauge. Define Ai{q) = 
J Ai exp{—iqx)Ax as the Fourier transform of Ai. On 
a circle of circumference i, configurations satisfying 
/ Aidx = Ai(0) = can be gauge transformed into con- 
figurations satisfying Ai{Q) = 21:111/ with n an inte- 
ger (the winding number). Therefore, when we consider 
a state that is an eigenstate of the quantized operator 
Ai(0), say, with eigenvalue 0, we should consider instead 
a linear combination of all states, each with eigenvalue 
27rn/e. The linear combination has to be gauge invariant, 
and the remaining arbitrary choice would be the phase 
between states with consecutive n. We call this relative 
phase 9 and call the vacuum of this phase the 0- vacuum: 
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Notice that if we write down the path integral from 
some state with winding n to some other state with wind- 
ing m by turning on Aiit), the winding number can be 
written as an integral: 



e 

2^ 



dxdt 



9Ai 



(5) 



here the limits of the integral denote the winding num- 
ber of the initial and final configuration. The vacuum- 
vacuum amplitude can thus be expressed as 
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here S in the exponent is just the ordinary action cor- 
responding to H and the scripts of the integral specifies 
the initial and final boundary conditions. The 0- vacuum 
description is thus equivalent to adding Ce to the La- 
grangian. 

Now we proceed to derive the physical consequence of 
the term. Consider a Maxwell Lagrangian with vacuum 
angle 9 at finite temperature 1//3. Taking Aq — 0, the 
Maxwell Hamiltonian is 
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dA,{q) 



dt 



E^i^^H.)p-E^. (7) 



Since Ai{q ^ 0) decouples from ^^(0) wc can calculate 
them independently. 9 only couples to the <; = sector 
as all operators at finite q have the same eigenvalue for 
states which differ by arbitrary winding. Let us focus on 
the partition function of the q ~ sector: 



Tre(e 



2ir 



d^ 



1^) (£|0 + 27rn) e'^™"")''; 



(8) 



the subcript 9 denotes that we only trace over the sector 
whose vacuum is the ^-vacuum, (p = eA{q = 0,t = 
0) is the initial configuration of the gauge field. Note 
that we have inserted 1 = ^ \£) where £ is the 

eigenvalue of {E^{q — 0)/eL) and \£) the eigenstate. The 
canonical conjugate pairs (x,p) can be determined from 
the Lagrangian with p ~ if we choose {eA^{q — 0)) 
as X it conjugates to {E^{q = 0)/eL). Therefore we have 

{(l) + 2mn\e) ^ cxp{i{(j) + 2TTm)e). (9) 

There is translational symmetry in m and n and the sum 
over m + n just gives an overall normalization constant. 
If we replace (m — n) by n, we have 

^^-"-f ^^e»(«+2-^)e-^^^ (10) 

If we sum over n first, we have 

Y^^^n(e+2.l) ^Y^g^e_^^^^^^ ^^^^ 
n m 



Physically, this means the effect of the uniform 9 term is 
to cause the average electric field to be quantized in inte- 
ger units of charges, but shifted by e9/2TT. This is a well- 
known result with open boundary conditions^, where one 
can imagine fractional charges at the end produce the 
electric field. With periodic boundary conditions it is 
less intuitive. 

If = TT, this would imply that the vaccum has two de- 
generate configurations characterized by J dxE ~ ±^e. 
The matrix elements between the two states become ex- 
ponentially small as L — >■ c», so we should think of this 
as a sponteneous symmetry breaking situation where the 
parity (P) and charge-conjugation (CC) symmetry are 
sponteneously broken by the electric field. The electric 
field would choose one direction and stay for a time pe- 
riod proportional to e^. 

In conclusion, there is indeed a real measurable differ- 
ence between Cg and C\/e, where with 6* = tt in the for- 
mer theory there will be huge ground state electric field 
at q ~ and the CC symmetry is sponteneously broken, 
whereas in the latter theory there will be no effect and 
the symmetry is preserved. 
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Now we turn our attention to three dimensions. We 
consider two settings without boundaries: the first is a 
3-torus, and the second is the 3-spherc. We restrain our- 
selves to consider only U(l) gauge fields. 

Abelian gauge field on a 3-torus: 

Since we need a periodic lattice to produce the topo- 
logical insulator, it is natural first to consider the world 
as a 3-torus. Again taking the gauge choice Aq = 0, the 
0-term in three spatial dimensions can be written as a 
difi^erence of the Chern-Simons term on the initial and 
the final states in the imaginary time direction: 

where i, j, k now run through only the spatial directions. 
One superficial difference to the situation in ID is that 
it seems all finite-g components contribute. However, as 
we require the initial and final states to differ from each 



other only by a gauge transformation, Afinai = ^initial + 
V0/e, we can see the integral on the three-dimensional 
boundary becomes a total derivative, 

j d?xd,{(j)djAke'^''/e)= j d?xd,{(l>B,). (13) 

Let us assume (j) only has a winding in the z direction, 
i.e., (j){x, y, Lz) — (t>{x, y, 0) = 2?™, then eq. ()13|) becomes 
27rn$ B where $b is the total fiux threading the torus in 
the z direction. Assuming $b = to$o with $o = hc/e = 
27r/e, we find 

/ d^xe^^'^^^df.A.dxA^ =. ^nm, (14) 

Js 6 

Thus, with m units of the fundamental fiux quantum in 
the z direction, the " 0- vacuum" consists of linear super- 
position of states with configurations satisfying J A^dz — 
27rn/c, where n is an integer. Since the Hamiltonian is 
still quadratic, we can calculate the relevant part of the 
partition function similar to the calculation in ID. The 
analog of eq. ([TU)) is 



V = LxLyLz is the world volume and we choose 
our conjugate variables to be (e j Azd'^x/LxLy) and 
(/ i^^d^x/eLz), with the eigenvalue of the latter labeled 
by i. We find that with a fix flux <^b = "m^Q in the z 
direction, the electric flux in the same direction is quan- 
tized: 

E^L^Ly ^ e{n - me/2TT) = ne - 6'$b/$^^ (16) 

with n an integer. 

Let us first take the strict T = limit. Here the ther- 
mal fluctuation of the magnetic flux is suppressed and we 
flnd that, the 0-term only has nontrivial effect if there is 
a flnite flux threading through. For 6 ~ tt, when we have 
an odd magnetic flux, the electric flux in the same direc- 
tion would be quantized in half units of e. The electric 
fleld goes to zero if the world volume goes to infinity, 
however. 

At finite T, the thermal fluctuation of the magnetic 
field can generate some finite fluxes, and we would have 
some effect even with _B = in average. For simplicity 
let us again set 6* = tt and consider = Ly = ^ 
L. If T ^ 1/L, The correlation function of the electric 
field would contain an extra term comparing to the usual 



Maxwell theory: 

{E{x)E{y)) ^ {E{x)E{y))\e=^ + (17) 

One can understand this constant correlation by imag- 
ining that half of the states in the ensemble have an 
odd number of magnetic fiuxes. The state with an odd 
number of fluxes would have a ground state electric field 
squared to (e/2i^)^, and the average is just a half of 
that. This extra part of the correlation function is long 
ranged, and can easily be distinguished from the Max- 
ell part. However, the magnitude again vanishes in the 
large L limit. Since it is not possible to have a 3D torus 
without embedding it in 4D space, these effects are of 
academic interests only. 

Before we end this subsection, we should note that 
from this calculation, it is clear that any local magnetic 
field will not produce any effect. Therefore, one would 
not see an electric field inside a solenoid, nor any charge 
at the end of it. 

Abelian gauge field on a 3- sphere and magnetic 
monopoles: 

Since we cannot have global nonzero magnetic flux in 
any direction in a 3-sphere, there will be no effect of the 
0-term. This is in contrast to the case with a magnetic 
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monopole, where it is predicted that there will be charge 
e{n — 0/27r) attached to it in a 6'-vacuum. This effect 
can be understood as follows: magnetic monopole is a 
singularity in terms of the abelian gauge field. Suppose 
we have a pair of monopole-antimonopole far away in 
a 3-sphere so that we have one fundamental flux going 
from one to the other. The geometry is now a 3-sphere 
with two punctures. From the calculation of the previous 
section we can see the electric flux threading from one 
hole to the other must be quantized, = e(n — 9/2-^), 
and we would attribute this as the charge of the magnetic 
monopole. 

Franz et. al. showed that there is Witten effect inside 
the topological insulatorj^ We emphasize here that this 
does not prove that a bulk 6'-term exists, as the Witten 
effect can also come from the " V^" theory, provided that 
we characterize it by ^? = inside the monopole. Given 
that a monopole can only live in a unit cell, and the band 
structure is absent in the unit cell, it is not unnatural to 
set ^ = inside a monopole. 

As a side note, if we consider nonabelian gague fields, 
the 0-term in general does have effect in a 3-sphere. This 
effect, however, is usually associated with the physics of 
instantons and is quite different from what we have dis- 
cussed. 



III. THE EFFECTIVE ^-TERM IN THE 
PRESENCE OF BAND ELECTRONS 

In this section, we investigate how the presence of 
fermions can alter the vacuum 0-angle. We review the 
topological band invariant which characterizes the topo- 
logical insulators in one and three spatial dimensions and 
show that they are related to the shift of 9. For the clear- 
ness of the formula, sometimes we take the units e = 1; 
i.e., we absorb e into the definition of the external elec- 
tromagnetic gauge fields. 

The topological insulator in ID is characterized by the 
polarizatioii^ 



P 



with 



dfc X ^ 



dfc 
27r 



tr (A) ; (18) 



{Unk\ 



\Un'k) 



(19) 



which is the so-called Berry's phase gauge field in mo- 
mentum space. \unk) is the periodic part of the Bloch 
wave function. 

P is forced to either take the value ^ or when there is 
charge conjugation (CC) symmetry. If the CC symmetry 
is preserved everywhere, then on the boundary there will 
be an odd number of zero modes. If the CC symmetry is 
locally broken in some way on the boundary, then there 
will be a 71 + 1/2 charge, where n depends on the detail 
of the local symmetry breaking. A cartoon showing the 
effect is depicted in Fig. ([T|). 
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(a) normal insulator 



C7) (+J l+O 



(b)topological insulator 



FIG. 1. Topological insulator in ID. The dashed line shows 
the edge. The envolope shown is the Wannier wave function 
of the electrons. 



Naively one might think this already shows that 9 is 
shifted by 2ttP: after all, the 0-term in ID is nothing 
but an energy term proportional to the electric field, in 
which the energy of dipoles, — J dxP • E, fits perfectly. 
However, we should note that normally this dipole energy 
arises from separating the charges to the boundary, away 
from each other. On a circle with uniform polarization, 
therefore, one would not anticipate such energy term is 
present. 

If we look back at how 9 change the physical property 
of the system, it comes in by adding a phase to the am- 
plitude between vacuua with different winding numbers. 
Specifically, 9 is precisely the additional phase of the am- 
plitude between vacuua with consecutive winding num- 
bers. In the presence of gapped fermions, this phase can 
come from integrating out the fermions, in other words, 
the dynamical phase the fermionic system obtains under 
a time-dependent back ground field. This phase has two 
contributions, one is just the time-dependent energy of 
the fermions and the other is the Berry's phase of the 
process. The phascshift from the energy depends on the 
time duration and is not just a function of the initial and 
final state; therefore it does not alter 9. Therefore, sim- 
ilar to the consideration in Refj^, we are led to consider 
the accumulated geometric phase of the band electrons, 
when the external field is slowly turned on. See Fig. HJa) 
for a cartoon of the procedure. 

First we shall consider how the single particle wave 
function change as we increase Ai uniformly. We have 



Unk 



{x)e 



ikx 



(20) 



which is the wave function in position basis, and u„fc(x) 
is periodic and satisfies 

((V - (fc + eAi)f + V{x)^ Ur,k{x) = EnkUnkix). (21) 

As we increase Ai uniformly to Ai -I- 77, the momentum 
k cannot change as it is fixed by the finite size L and 
the periodic boundary condition. On the other hand, 
following cq. (HI]), Unk{x) changes as 



-lnk{Ai +7?) = Un{k-e7j){Al), 



(22) 



which is just a corresponding shift of k by —erj. if 
erj = 27r/L, the system returns to its original state, but 
in a different gauge (i.e., with winding number different 
by one.) Notice that while Unkix) goes to the next ava- 
iable value on the left, the k in the exponential stays the 
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(a)thc procedure 



(b)spcctral flow 



FIG. 2. (a) A flux is slowly threaded through. <E> = and 
$ = 27r are the same physical state related by a gauge trans- 
form. We calculate the Berry's phase of the process, (b) 
During the process, at every allowed momentum by the peri- 
odic boundary conditions, the energy and the periodic part of 
the wave function moves slowly to the values of the state to 
the left, according to eq. (|22|l . When a full flux is threaded, 
each one of them takes the eigenvalue and the eigenvecotor 
of the one at its left. Note that the momentum quantum 
number k, however, does not change. When we sum over the 
Berry's phase contribution from all the single particle states, 
it becomes an integral over the entire Brillouin zone. 



same. The electronic wave function is therefore different 
from its starting state. Nevertheless, as discussed fur- 
ther below, if we include the gauge field, the final state 
differs from the initial state by a large gauge transforma- 
tion, and the Berry's phase accumulated in the process 
is exactly what we want to calculate. 

As a side note, the situation is similar if we put elec- 
trons on a lattice which couples to the gauge field via 
Periels substitution. The single particle eigenfunction 
can be written as Tpnk = J2i '^nk,m e'xp{ikxi) |m, a;,), with 
Unk now a vector in the orbital space. With an increase 
of Ai, only Unk changes. 

Now we are ready to calculate the accumulated Berry's 
phase of the band electrons under the process, where the 
winding of the gauge field is increased by one: 

/.27r(«+l)/e ^ Q 
/ dAl(0)(*ol^ 



aAi(o) 



' 27rn/e 

r27r(«-|-l)/e ^ Q 
I I dAi(O) 2J (V'k.al ^~ l^kia) 



2-Kn/ e 



aAi(o) 



l-k,+27T/L Q 
^ j d/c (Ufecl — 



/" dfc {Uak \ ^ \Uak) 

Jbz (^^^ 



aGocc 



= -27rF. 



(23) 



In the second equality, we wrote the derivative acting 
on the Slater determinent as a sum of derivatives acting 
on single particle wave functions. In the third equality 



we then plug in the dependence of the wave functions, 
and change variables to k. Whenever Ai{0) increases by 
27r/e, each Unk reaches the next allowed eigenstate to the 
left by the periodic boundary condition (without actually 
changing the momentum eigenvalue.) As we sum over all 
the integral of eigenstates at different allowed fc's, the 
whole Brillouin zone (BZ) is covered exactly once and we 
reach the fourth equality. 

While we calculate the Berry's phase for process where 
the winding number of the initial and final states differs 
by one, evidently the phase is proportional to the dif- 
ference of the winding number in general. Therefore, it 
leads to a shift of 9 by (/)Borry Specifically, if the vacuum 
has an vacuum angle 9 = 0, we see that the topological 
insulator in ID will be described hy 9 ~ 

A few remarks are in order. Firstly, this calculation 
is good for a finite-size system, where both the lattice 
spacing and the length of space are finite. Despite that 
the eigenstates in such a system would be discrete points 
in the BZ, the whole BZ is covered by the integral and 
there is no finite-size effect. Secondly, one might notice 
that the fermionic wave functions of the initial and final 
state are different, and it seems that our process does 
not form a close loop as usually required by a physical 
(gauge-invariant) Berry's phase. This does not invalidate 
our calculation, however, since the initial state and the fi- 
nal state are related by a large gauge transform. Once the 
convention of the phase of the initial state is determined, 
the phase convention of the final state is also determined 
via gauge-transforming the initial state to have the de- 
sired winding number. The residual gauge degree's of 
freedom without altering the external gauge field (which 
is just an arbitrary phase of the final state) can be ab- 
sorbed into the vacuum angle 9 in the absence of the 
fermions. 

Now we turn to the 3D strong topological insula- 
tor. The 3D strong topological insulator, defined un- 
der TR symmetry, is characterized by the band structure 
invariant^ 



^ J d3A:e'^''^Tr( - ^^ A A A) = n. 



(24) 



Similar to the case in ID, a 0-term is generated if we 
can find a Berry's phase of the electrons which is pro- 
portional to the difference of the winding number of the 
initial and final gauge configuration. More specifically, 
the procedure is as follows: we apply a constant finite 
magnetic field, say, in the z direction on the 3-torus. We 
then slowly change the gauge field in the same direction 
uniformally until the final state is connected to the ini- 
tial state by a large gauge transform. Then we apply 
the magnetic field in some other direction and repeat 
the calculation. We can also consider procedures such as 
applying magnetic field in z direction and changing the 
gauge y direction; the phase of this process leads to a 
term oc EyB^ in the effective Lagrangian. In general, we 
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therefore expect the full effective theory to take the form 
/:cff = ^a,j£;,Bj. (25) 

The 0-terni is the isotropic part of the effective La- 
graigian: 



47r2 
3? 



(26) 



Two things are different from our calculation in ID: 
firstly, the state at a given kz is already a many-body 
wave function. Secondly, we have to calculate everything 
in a finite (but maybe small) magnetic field. In principle 
one can go over all the Laudau levels at a given and 
sum up their Berry's phases, but in practice this is not 
easy to do as the wave functions are not perturbative 
in B. It turns out that the density matrix perturbation 
theory introduced in Ref.— is suitable for this calculation 
with an extra trick as will be described below. 

Before we dig into the calculation, let us clarify that 
our calculation, despite taking advantage of the same for- 
malism, is distinct from Refji^,. There they first calcu- 
late the current flowing through the bulk as they vary the 
Hamiltonian with time under a small magnetic field, then 
they relate the time integral of the current to the polar- 
ization. While the uniform 0-term in 3D with boundaries 
would give rise to a magneto-polarization effect, the con- 
verse cannot be said. As mentioned in the introduction, 
both Cg and £ve can produce this effect, so a deriva- 
tion of the effect does not distinguish between the two 
theories. To our best knowledge the following calcula- 
tion is the first demonstrating that it is indeed Cg which 
describes the topological insulator. 

Suppose we apply a small magnetic field along the z- 
direction, B = Bz. Following the calculation in ID, we 
calculate the Berry's phase of the process: 



^Berry 



^/dkz(*Bl|;l*B); 



(27) 



here ^I*^ denotes the slater determinant of the 2D elec- 
tron wave functions, for a given kz in the magnetic field. 
Analogous to the case in ID, the derivative is understood 
to be taken only on the periodic part of the Bloch wave 
function. 

Despite that the integrand can be written as a sum 
over single-particle wave functions, we immediately no- 
tice that it cannot be expressed as a function of the single- 
particle density matrix. This is due to the fact that the 
wave function depends on the vector potential which is 
gauge dependent. One easy way to realize the fact is to 
consider a change of phase in the wave functions. The 
integrand is not invariant (the integral as a whole, on the 
other hand, is invariant modulo 27r) whereas the density 
matrix remain unaltered under the transformation. 

We can, however, express the integral as a whole in 
terms of the density matrix by the following trick. The 



accumulated Berry's phase can be expressed as an inte- 
gral of the Berry's curvature in one extra dimension using 
Stokes theorem, with the region of integration bounds 
by the original kz integral. The Berry's curvature can 
now readily be expressed in terms of the density matrix 
extended into the extra dimension, which is chosen con- 
tinuously but otherwise aribitrarily with the constraint 
such that on the boundary, we have the original density 
matrix. We therefore have 



dkz(*Bl^|*B) 

as c'^z 



*B 



d2ke"^Tr(p9„p9^p); 



(28) 



^ B 



is the 2D electron many-body wave function, 
which in addition to being a function of fc^, has been 
extended to some extra direction kw p = X^i V'i^ i^i 
is the extended 2D single particle density matrix, where 
ipi ) is a 2D single particle eigenstate in magnetic field 



B. a, 13 run through two directions which is spanned by 
kz and A:^,. The trace sums over both the band and the 
(x, ?/)-position basis. As mentioned above, the integral is 
chosen to be performed on the area such that the bound- 
ary is at {kz, kw = 0) and the density matrix has known 
values. It is straight forward to show the second equality, 
and the derivation is provided in the appendix. Different 
choices of density matrices inside the boundary can only 
alter the integral by multiples of 27ri. To avoid cluttering 
of the equations, in the following we omit the tilde for 
the extended objects when there is no ambiguity. 

Then following the formalism in RefJ^, we take the 
large size limit and expand p to linear order in B. As 
discussed there, the density matrix in real space basis 
can be decomposed into two parts, one of which is trans- 
lationally invariant: 



Pr 



/9ri,r2 exp(-ii? • (fi X r2)/2). 



(29) 



where /Cri,r2 denotes the density matrix in position ba- 
sis, and p is translationally invariant. While the other 
part seems to be affected by the infinite range of r, in 
our expression three p's appear together and the combi- 
nation is short-ranged and can be expanded in B. It is 
thus straight forward to expand p explicitly and calcu- 
late. (Please see appendix for details of the calculation.) 
Up to first order in B, the result is 



0Borry 



d2fc' 



(fk 

s Jbz (27r 



-idaTr^dppok'il - Pok')pk'Pok' ~ h.c.) ;(30) 



a,/? span kz,kyj and 7, i5 span k^,ky. The integral of 
k' is performed on the 2D Brillouin zone in xy-plane. 
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Pk' = (^'IpI^') is the translationally invariant part of 
the density matrix at a given [k,kw) and pofc is the den- 
sity matrix in zero field. J^fm is the nonabehan Berry 
cmvatme of the occupied bands, 



J~ dfi-^i' 9i/jA^ i[jA^^ jAi/^. 



(31) 



Notice that in eqn. ([30| . the tensor structure in the 
first and the second term is different and we can rewrite 
the first term using the total-antisymmetric tensor in 4 
dimensions: 

0Borry = 4>I + 4'A 

'^^ = ^ / d4ke'^^^'^Tr(J-,faJ-ed) (32) 



1^ 



(33) 



SxBZ 



dpk' 



Ma0 ^ Tr ydaPok'i^ - Pok')-g^Pok' - h.c.j . (34) 

We have explictly expanded the second term to first order 
in B. a, b, c, d runs through all directions. Both integrals 
are total derivatives and we can integrate back to the 
boundary which is the original 3D Brillouin zone: 



B 



87r2 
47r2 



d-^ke'^'^TiiAadbAc - i=-AaAbAc) (35) 



(36) 



(/)/ is isotropic, i.e., independent of the direction of the ap- 
plied magnetic field. 0^, on the other hand, is anisotropic 
in the sense that if we do the same calculation for the 
magnetic field in x or y direction, the result in general 
would be different. Now we consider the gradual gauge 
transform in the i-direction and the magnetic field in the 
j-direction, the same calculation still goes through, pro- 
vided that we take a, /3 in the i-direction and the extra 
direction, and 7,(5 in the directions prependicular to the 
magnetic field. We get 



47r2 



(TkNh. 



(37) 



In terms of the effective theory, this means that the ef- 
fective Lagrangian not only contains E ■ B, in general we 
have J2ij OitjE^Bj, where 



d^fc 
(2^ 



e-bc -Lr^r (AadbAc - l\AaAbA,)h, + M,. 



(38) 

By calculating the Berry's phase of these processes, not 
only do we get the coefficient of the topological term but 
we also get a part which is a physical response which 
agrees with Ref.— . In general Ma also contributes 
to 0. If TR symmetry is present then Mj^ = and we 
see that the vacuum angle is shifted by tt in the pres- 
ence of the strong topological insulator. We stress again 



that the calculation present here shows directly that the 
vacuum angle B is shifted in the presence of the elec- 
tronic band structure whereas the previous calculations 
only show that one can get current responses when one 
smoothly varies the Hamiltonian. Physically, our result 
predicts that there will be a half-charge electric flux if we 
put a strong topological insulator on a 3-torus with an 
odd number of magnetic flux, as described in the previ- 
ous section; whereas in previous derivations, it is unclear 
if one can observe anything without either a boundary or 
a change of the Hamiltonian. 



IV. SUMMARY 

In this paper, we flrst show that there is a measur- 
able difference for an effective theory containing either 
or iZvei in one and three spatial dimensions without 
a boundary. Specifically, with /Ig, in ID there will be an 
electric field 0e/27r in the ground state. When 9 = tt the 
electric field can be in either direction and the CC sym- 
metry and parity symmetry is spontaneously broken. In 
3D, the same effect can be found on a 3-torus, but since 
it is the electric flux which is proportional to 6e in 3D, 
this effect vanishes in the thcrmaldynamic limit. 

We then go on to show that the topological insulators 
in ID and 3D can be described by Cg instead of C-^e. 
While the expression of agrees with previous results, 
to our knowledge this is the first derivation which distin- 
guishes between Cg and C\jg. 

We thank N. Nagaosa, X. G. Wen and F. Wilczek for 
helpful discussions. We acknowledge the support of NSF 
under grant DMR 0804040. 



V. APPENDIX 

A. Path integral formulation for the 0-term 

In the main text, we derive the physical consequence of 
the ^-term using the notion of 0- vacuum, which is similar 
to a Hamiltonian formalism. One may wonder why we do 
not directly carry out the path integral. The first reason 
is that the quantization is not obvious if we just do the 
euclidean path integral as done below. The second reason 
is that if we calculate the fluctuation of the electric held 
at flnite temperature, a naive calculation would give us 
a sum of negative values, which does not make sense. 
It turns out that for a free theory the position space 
path integral is ill-behaved and a positive finite term is 
expanded as an infinite negative sum. A similar situation 
occurs when one calculates the ground state energy of the 
bosonic string using mode expansion. Let we start right 
from the Lagrangian 



1 

4^ 



ee 



(39) 



using the gauge — 0, for the q — Q sector at finite 
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temperature 1/ (3 we have the partition function 

1 27rn 2 

(40) 

where cj^ = 2-Ki/ 13 are the Matsubara frequencies. Again 
the finite frequency part decouples and the zero frequency 
part agrees with eqn. ljlOp if we integrate I first instead: 




Nevertheless, it is hard to see from this form that 6 cor- 
responds to a quantization condition for the electric field. 
Without reversing the i integral, another way to see the 
0-dcpendcnce is to calculate the expectation value and 
the fluctuation of the electric field. For 9 = tt, the expec- 
tation value would vanish and we can only rely on the 
fluctuation. 

If we calculate {\E^ (0) P) , however, we would encounter 
a problem in the path integral as now all flnite frequency 
part contributes and their sum seems to be inflnitely neg- 
ative: 




f2TTn\' 



If we compare this to what we would have got using the 
Hamiltonian formalism, 



It seems we have to require 



27rny 1 



(43) 



(44) 



for the two expressions to agree. We can understand 
this equality by thinking of the left hand side as the zeta 
function at zero, C(0), written in a series. While the series 
is divergent at zero, the zeta function is well-defined and 



is indeed — i. 



The function J^n ^^^n is related to the elliptic 0- 
function. If one subtracts the fiuctuation at = from 
the expression and calculate at /3 — oo, one recovers that 



(|i?i(0)p)-(|#(0)p)|,.o = iV - 



2n 



(45) 



which implies the quantization. 



B. Derivation of eq. ([28]) and eq. ([30]) 

Here we show explicitly the derivation of the second 



equality in eq. 



First we note l^fs) = det(?/'i)j a-nd 



Now we plug in p = ^ - l^i) to the right hand side 
of the second equality, we have 

Tr (pdapdpp) 

= E 1^') ((^" i'^^)) (^^1 + 1^*) (^^d) 

ijk 

+Y.{d^midp\^..)y, (47) 

i 

in the derivation we have taken advantage of the fact that 
(V'jIV'j) = Sij and thus (a„ {■i/)j|) l^/'j) = - {tpi\ IV'j)- 

Contract both eq. and eq. (|T7)) with e"^, we can 
see that they agree. 

In the following we apply the magnetic field in the z- 
direction and take p as a function of kz , and r lies in the 
xy-plane. We take fi, = e = 1. p has the following matrix 
elements up to first order in B: 



('0mfc|p|V'm'fc') 
{■>Pnk\p\'>Prak') 



^kk' 



jS.nn' 



■^Skk'Be^^TjS,mm' 

. -yS (V'nfcl {d^POk, dsHk} \lpvik) 

''^'V2 " Enk 



Err 



{ljjnk\H'f, \lprnk) 



(48) 



Enk — E,nk 

note that the momentum here is two-dimensional and ev- 
erything has implicit k^, k^ dependence, n, n' are indices 
for occupied bands and m, m' are for empty bands. T is 
the nonabelian field strength for the Berry's phase gauge 
field defined from the unoccupied bands: 

, d 



(49) 



In the following computation one would find these ex- 
pressions useful: 



-F Ul^.l 



T 

■J LLV. 



-^E (V'nfcl C^/x \lpTnk) {'ipnikld^ \i>n'k) 
m 

-i E (^'mfel <9^ \ Aik) {i^nk\d„ IV'm'fc) 



Xij) 



-{p^v)\ 
e^''^"Tr(J-^,,J-,A) 



(50) 



Note that in the expression for the Berry's curvature 
we use the whole Bloch wave function jV') instead of the 
periodic part \u) but here it makes no difference. 

Now we start from eon. ([25)1 . Write explicitly in posi- 
tion basis, we have 



=E^"(V^^I^/3l^')- (46) Tr((a„p) p (dpp)) ^ J dndr,dr,id^p,,)p,,idppsi) 
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exp (^2 - ri)^{r3 - ri)sj 

(51) 

where in the second equahty we taylor-expand in B, keep 
up to first order and go back to momentum space. We 
have also taken the infinite-size hmit and make the sum 
of discrete momenta an integral. The trace on the right 
hand side traces over only the band indices. 

The remaining task would be to plug in p and calcu- 
late explicitly to first order in B. One thing to notice is 
that when taking derivatives of p, it acts not only on the 
matrix element but also on the basis. It is also useful to 
note that duPo only has non- vanishing matrix elements 
between the original occupied and empty states. 

As we can see from ean. (j48p . the inter-gap and intra- 
gap matrix element of p look pretty different. Let us 
denote the former as p' . p' contributes only through the 
first term in the right hand side of ean.([5T|): since p' is 
already first order in B the remaining p can be replaced 
by pq. Let p' = A+A'' with A = (1 — po)p'po (that is, A is 
the matrix clement connecting occupied bands to empty 
bands and vice versa for A^), after explicit calculation, 



similar to eq. (|47)) . we have 

e"^ {Tiida^p'podppo) + Tiid,,pop'dppo) + Tr{d^popodpp')) 
= -da i{n\ d/3 \m) A,„„ - c.c.) 

= 9„Tr {dppoil - po)p'po + h.c.) ; (52) 

|n) is the short hand notation of \ipnk) and repeated in- 
dices are summed over. 

Now that the inter-gap matrix elements are dealt with, 
the remaining part of the first term can also be expanded 
and calculated: 

e"^Tr(9,pp9^p)|remai„i„g = e"^(^Tr(J-„^) 

+ 0{B^)). (53) 

The first term on the right hand side is proportional to 
B^ and is similar to the polarization in ID. 

The only remaining part is the second term in eqn. (|5ip . 
This part proves to be somewhat tricky to calculate as 
one has to manually group terms into expressions of J- 
and J^. Nevertheless, it is otherwise straight forward and 
one gets 

e'^'^e^^Trid^daPoPodsdppa) = 

Tr ^-TapJ'jS — -^J-apJ^'jS + ^Sa^-fP^ ■ (54) 

Combining egn.dS^. can. ((5^ . and cqn. ([M)) . and with the 
help of eqn. ([50]) we get egn. (PO)) . 
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